arXiv:1508.03941vl [cond-mat.str-el] 17 Aug 2015 


The antiferromagnetic cross-coupled spin ladder: quantum fidelity and tensor networks approach 
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We investigate the phase diagram of the cross-coupled Heisenberg spin ladder with antiferromagnetic cou¬ 
plings. For this model there have been conflicting results for the existence of the columnar dimer phase, which 
was predicted on the basis of weak coupling field theory renormalisation group arguments. The numerical 
work on this model has been based on various approaches, including exact diagonalization, series expansions 
and density-matrix renormalization group calculations. Using the recently developed tensor network states and 
ground-state fidelity approach for quantum spin ladders we find no evidence for the existence of the columnar 
dimer phase. We also provide an argument based on the symmetry of the Hamiltonian which suggests that the 
phase diagram for antiferromagnetic couplings consists of a single line separating the rung-singlet and Haldane 
phases. 

PACS numbers: 75.10.Jm 


I. INTRODUCTION 

Geometric frustration induced by competing interactions in 
quantum systems produces increased quantum fluctuations re¬ 
sulting in a range of distinct quantum phases with potentially 
rich phase diagramSfi Such systems can be low-dimensional if 
there are a sufficient number of competing interactions, which 
need to be comparable in magnitude to induce frustration. The 
example which we examine here is the cross-coupled two- 
leg Heisenberg spin ladder. Spin ladder systems are of con¬ 
siderable importance in condensed matter physics and have 
thus been widely investigated, with a number of known ladder 
compounds realised in experiments^ The cross-coupled spin 
ladder has Hamiltonian 

H(J\\,J±,Jx) - -fll l,i • S i,,'+i -H S2,i ■ S2,i+l) 

i 

+ Jl. X ^ 1.' ■ ^ 2,1 

i 

+ Jx ^ 2,1+1 + ^ 2,i ■ s l.i+i), (1) 

i 

where S i,, and 82,1 are spin-1/2 operators on rung i and legs 
1 and 2. The parameters yy, J± and Jx are the intra leg, in¬ 
ter leg (rung) and cross exchange couplings, respectively. The 
spin ladder couplings are depicted in Fig. [T] We are most in¬ 
terested in competing antiferromagnetic couplings, for which 
there have been conflicting results for the precise nature of the 
phase diagram. 

The cross-coupled ladder model in Eq. ([TJ has been in- 
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FIG. 1: The cross-coupled two leg spin ladder. 


vestigated, using numerous methods, for over already two 
decades4^ After some debate it now appears settled that the 
phase diagram for antiferromagnetic couplings consists of a 
phase boundary separating rung-singlet and Haldane phases. 
As emphasized^S the entire region above the transition line is 
in the rung-singlet phase, being continuously related to the ex¬ 
actly solvable line Jx - J\\ along which an exact rung-singlet 
state is the ground-state in this phase. The entire region below 
the transition line is in the Haldane phase, being continuously 
related to the point iJxlJ\\ - ^,Ja.IJ\\ - 0) at which the ladder 
is equivalent to a spin-1 chaini^— 

The debate centred over the existence or not of a colum¬ 
nar dimer (CD) phase intermediate between the rung-singlet 
and Haldane phases which was predicted on the basis of weak 
coupling field theory renormalisation group arguments^i This 
phase is between J^^ 2Jx - {5Jl)/{7rJ\\) and 7^ - ^Jx - 
{J^)l{7T^J\\)ii^ Subsequent numerical work motivated by this 
prediction disputed the existence of the CD phase On 
the other hand, further numerical work provided evidence 
supporting the existence of the CD phase,2ii2^ In particular, 
the CD phase was claimed to exist in the range 0.373 < 
J± < 0.386 for Jx = 0.2^ This point was further ad¬ 
dressed by more extensive density-matrix renormalization 
group (DMRG) calculations in the range 0.36 < Tj. < 0.4 
for Jx = 0.2, from which it was concluded that the CD phase 
does not appear^ The most recent work highlighted the fun¬ 
damental limitations of the DMRG method when applied to 
the cross-coupled ladderi^ Nevertheless, the existence of the 
CD phase in this model is beyond dispute for ferromagnetic 
couplings, as confirmed by a number of studiesi22i22i2^ 

An efficient tensor network algorithm has been developed 
for spin ladder systems which generates ground-state wave 
functions for infinite length quantum spin ladders.— In gen¬ 
eral, tensor networks have proven to be a convenient means 
to represent wave functions of quantum many-body lattice 
systemsr^ and in particular, to calculate quantum fidelity as 
a measure of quantum state distinguishability. The ground- 
state fidelity per lattice site based on tensor networks has thus 
been established as a powerful tool to study quantum phase 
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transitions regardless of what type of internal order exists in 
the quantum states. It has been demonstrated to directly (i) 
capture phase transition points due to changes of ground-state 
wave functions in the vicinity of quantum critical points22r— 
and (ii) detect degenerate ground-states due to spontaneous 
symmetry breaking.—— The quantum fidelity approach can 
thus be utilized to provide evidence for the existence or not 
of the CD phase, because by its definitionri22i2i2^ the CD 
phase is to be induced by a spontaneous breaking of lattice 
translational symmetry. Motivated by the conflicting results 
for the existence of the CD phase in the antiferromagnetic 
regime and the advantages of the quantum fidelity and tensor 
networks approach, we have examined the cross-coupled lad¬ 
der model, paying particular attention to the previously stud¬ 
ied line Jx - 0.2. 

The paper is organised as follows. In Sec. II we outline the 
numerical procedure for the tensor network based quantum 
fidelity approach for the cross-coupled spin ladder and define 
the relevant order parameters. Our numerical results using this 
method are also presented. Concluding remarks are given in 
Sec. III. 


II. NUMERICAL PROCEDURE AND RESULTS 

The ground-state fidelity and tensor network states ap¬ 
proach has been developed to study the criticality and phase 
diagrams of ladder models in Ref. We refer to that paper 
for technical details. The tensor network representation suit¬ 
able to describe the ground-state wave function of an infinite- 
size spin ladder is a particular case of projected entangled-pair 
states (PEPS)^ adapted to the geometry of the spin ladder sys¬ 
tem. For the two-leg ladder one needs to update four different 
four-index tensors. In this prescription the updating procedure 
is closely connected to the infinite PEPS algorithm^ and the 
translationally invariant MPS algorithm!^ The computational 
cost scales as D® where D is the truncation dimension. For the 
two-leg Heisenberg spin ladder it was demonstrated that the 
tensor network approach could very accurately determine the 
ground-state wave function, with truncation dimension D = 6 
outperforming the DMRG results— 


A. Ground-state fidelity per lattice site 

1. Phase transition point 

For the cross-coupled spin ladder in Eq. ([1]) we set 7|| = 1, 
fix Jx and vary J± as control parameter. With a randomly 
chosen initial state in the tensor network representation, the 
updating algorithm of the tensors until the saturation of the 
energy in Ref. allows us to obtain the ground-states of the 
spin ladder. For a given with Jx - 0.2, we calculate the 
ground-state \ift(J±)). In general, one can consider the quan¬ 
tum fidelity F{\i//(J±)}, |0)) = KiA(./±)|0)| which is the overlap 
function between the calculated ground-state \il/{J±)} and an 
arbitrary reference state |0). The quantum fidelity scales as 
F(\ip(J±)), \4>)) ~ d(\iJ/(J±)), where L is the system size. 



FIG. 2: (color online) The ground-state fidelity per lattice site 
d(Ji_, J'j) as a function of and J'^ for the cross-coupled spin lad¬ 
der with coupling Jx = 0.2. (a) A continuous phase transition point 
0.383 is identified as a pinch point J\) on the fidelity sur¬ 
face. Two different phases are thus identified: the Haldane phase for 
J^ < Jl and the rung singlet phase for Jj^ > J^. (b) The correspond¬ 
ing planar view. 


Following Ref. [HI the fidelity per lattice site can be defined 
as 


lnd{tf/{J_L)},\4>}) = lim 

L—»oo 


lnF(|iA(7^)),|0» 

L 


( 2 ) 


It satisfies properties inherited from the characteristics of the 
fidelity F{\\l/{Jj_)), |0)), namely (i) normalization r/(|i/'(/±)) = 
10), |0>) = 1 and (ii) range 0 < ^f(|0(/±)), |0» < 1. 

When the arbitrary reference is chosen as a ground-state of 
the spin ladder, the ground-state fidelity per lattice site has 
been demonstrated to be a universal marker for detecting the 
existence of quantum phase transitions, which are character¬ 
ized by a pinch point on the fidelity surface, irrespective of the 
dimensionality of the system—^ For our system, then, one 
can choose the reference state |0) = |0(7j^)) to find a pinch 
point, i.e., a phase transition point. In this case, the fidelity per 
lattice site satisfies one more property, namely (iii) symmetry 
d(JjL, J'i) - d{J'_^, J±). In Fig. [2] we plot the fidelity surface 







































3 


d(\{l/(JAy), as a function of the rung coupling com¬ 

puted with truncation dimension D = 24. The fidelity surface 
satisfies the symmetry property for exchanging 7^ and jy We 
find a pinch point at (7 _l, 7^) = (0.383,0.383) in the param¬ 
eter range 0.25 < J± < 0.5 with 7x = 0.2 and 7|| = 1. The 
pinch point at (7j., 7^) = (0.383,0.383) implies that a contin¬ 
uous phase transition occurs at 7^ ^ 0.383. Most importantly, 
only one pinch point in the fidelity surface is observed in the 
parameter range, which indicates that there are two character¬ 
istic phases. However, the CD phase was predicted to be in¬ 
termediate between the rung-singlet and Haldane phasesjii In 
other words, the existence of the CD phase implies that there 
should be two pinch points in the parameter range. Conse¬ 
quently, the observed single pinch point on the fidelity surface 
indicates the non-existence of the CD phase in the parameter 
range. 


2. Detecting non-degenerate ground-state 

One can make a further check on the existence of the CD 
phase, which has not been applied to ladder models before. 
According to the Landau theory for quantum phase transi¬ 
tions, if a phase can be characterized by a local order param¬ 
eter, the phase is induced by a spontaneous symmetry break¬ 
ing and, as a consequence, the system should have degener¬ 
ate ground-states in the phase. If this theory is applicable for 
the CD phase, the CD phase should have degenerate ground- 
states. Because the CD local order parameter— is given by 

OcD “ 9 ^ ^ a,i+\ “ S a,i+\ ' S a,i+l) , (3) 

a=l,2 

the CD phase should be responsible for two degenerate 
ground-states. It vanishes in the rung-singlet and Haldane 
phases. This local order parameter should be nonzero in the 
CD phase and is thus a measure of the existence of the CD 
phase. 

To determine the existence of the CD phase, then, we 
need to detect degenerate ground-states. The ground-state fi¬ 
delity per lattice site has been utilized in detecting degenerate 
ground-states for spontaneous symmetry breaking—— To ap¬ 
ply this scheme here, for a fixed system parameter, we need 
to calculate ground-states repeatedly with n differently cho¬ 
sen initial states for the tensor network representation. Par¬ 
ticularly, the system parameters are chosen to be 7j. = 0.379 
with 7|| = 1 for 7x = 0.2, because these values were sug¬ 
gested to be in the CD phase— For our numerical study, the 
initial states have been chosen randomly and used to calculate 
ground-states with the truncation dimension D = 12. The ref¬ 
erence state \(p) has been also chosen randomly and has been 
kept to calculate the fidelity, which allows us to distinguish 
different ground-states with the same saturated energy in the 
numerical calculation. In Fig. [3ta) we plot the fidelity per lat¬ 
tice site d for the number of initial trial states n. Twenty initial 
trial states are here presented to produce the same ground-state 
because there is only one observed value for the fidelity per 
lattice site d. As a result, for this set of system parameters, the 
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FIG. 3: (color online) Fidelity per lattice site d between the ground- 
state and a reference state. The horizontal axis denotes the number 
n of the random initial state used to construct the ground-state. Also 
show is the corresponding local order parameter Ocd- (a) Antifer¬ 
romagnetic parameter values 7|| = 1, 7x = 0.379 and 7x = 0.2. In 
this case both the fidelity per site and order parameter exhibit one 
distinct value indicative of a non-degenerate ground-state, (b) Ferro¬ 
magnetic parameter values 7|| = 1, 7x = -1.6 and 7x = -0.7. This 
point is known to be in the CD phase for ferromagnetic couplings. 
In the CD phase it is clearly observed that there exist two degenerate 
ground-states corresponding to the two values of d and Ocd- 


system does not undergo any spontaneous symmetry breaking 
and has a non-degenerate ground-state, which implies non¬ 
existence of the CD phase. The vanishing local order param¬ 
eter Ocd in Fig. [3ja) also manifests the non-existence of the 
CD phase. 

In contrast to the case of antiferromagnetic couplings, the 
CD phase is known to exist for ferromagnetic couplings (7x 
and 7x < 0)— In order to compare with our results for 
antiferromagnetic couplings, we calculate degenerate ground- 
states for the parameter values 7|| = 1, 7x = -1-6 and 
7x = -0.7. These values are known to be deep in the CD 
phase for ferromagnetic couplings— In Fig.[3jb), we plot 
the fidelity per lattice site d for the number of initial trial states 
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FIG. 4: (color online) The string order parameters Oan and Oeven 
and the local order parameter Ocd as a function of for the cross- 
coupled spin ladder with coupling Tx = 0.2. Thin lines are guides 
to the eye. The system clearly exhibits the Haldane phase for 7^ < 
0.383 and the rung singlet phase for Jx_ > 0.383. The local order 
parameter Ocd is zero, indicating the absence of the CD phase. 

n. In the CD phase it is clearly observed in the same way 
that there exist two degenerate ground-states corresponding 
to the two values of d. In fact, for a large enough number of 
random initial state trials, the probability P{n) of each degen¬ 
erate ground-state approaches 1/2, i.e., Y\m„^ooP{n) - 1/2. 
This gives an explicit detection of the two degenerate ground- 
states for the CD phase which originate from a Z 2 sponta¬ 
neous symmetry breaking. Moreover, as expected, the local 
order parameter Ocd exhibits two distinct values that have the 
same amplitude and opposite sign, i.e., where 

the local order parameters and are calculated from 
the two degenerate ground-states |i/'g') and |iA®), respectively. 
The relation between the two local order parameters mani¬ 
fests a Z 2 symmetry breaking for the two degenerate ground- 
states. Consequently, the sharp contrast between Fig.|3a) and 
Fig. |3lb) lends further weight to there being no CD phase for 
antiferromagnetic couplings. 

B. String order parameters 


In the previous subsection, we obtained the ground-states 
as a function of for Ty = 1 and 7x = 0.2 with trun¬ 
cation dimension D = 24 to obtain the ground-state fidelity 
per lattice site in Fig. |2] In Fig. H) we calculate and plot the 
order parameters from the ground-states in Fig. |2] Note that 
the local order parameter Ocd is zero for the parameter range. 
Also, both the even/odd string order parameters Oeven!Oodd 
vanish at Jl ^ 0.383. The string order parameters clearly 
identify the Haldane phase for < 0.383 and the rung singlet 
phase for > 0.383. As is detected from the pinch point on 
the fidelity surface in Fig.|2l the behavior of order parameters 
shows that a continuous phase transition occurs at 7^ ^ 0.383. 
Consequently, the order parameters in Fig. |4]show further that 
there is no evidence for the CD phase. 

III. CONCLUSION 

We have examined the phase diagram of the cross-coupled 
Heisenberg spin ladder in Eq.([T]i with antiferromagnetic cou¬ 
plings using the recently developed tensor network states and 
ground-state fidelity approach for quantum spin ladders.— 
This method has been demonstrated to be particularly effec¬ 
tive for identifying quantum phases and transition points. We 
examined the particular coupling region 0.36 < J± < 0.4 for 
Jx = 0.2 for which there have been conflicting results for the 
existence of the CD phase. We identified a continuous phase 
transition at the value Jf ^ 0.383 and clearly observed the 
Haldane and rung-singlet phases. However, we did not find 
any evidence for the existence of the CD phase. This conclu¬ 
sion was further supported by the use of a related spontaneous 
symmetry breaking method, which we tested in the vicinity of 
the critical point and for comparison in a region with ferro¬ 
magnetic couplings where the CD phase is known to exist. 

A particular symmetry of the Hamiltonian in Eq. ([T]i poses 
a strong constraint on the nature of the phase diagram for this 
model. It was observed that H{J\\,J^,Jx) - H{Jx,Ji_,J\\)i^ 
which is readily seen to be evident by inspection of the Hamil¬ 
tonian. Defining the variables yj_ - J±/J\\ and jx = Jx/J\\ 
this symmetry implies the “duality” relation //(l,yj_,yx) = 
yxHil,y±lyx, 1/yx). Thus the phase diagram for yx < 1 maps 
to the valuesi2£ 


Order parameters for distinguishing between the rung- 
singlet and Haldane phases are nonlocal string order 
parameters— Eor two-leg spin ladders, we investigate the 
string order parameters Dodd and Devon defined byiSil 


<5p,exp 


ITT 


./ ^ 

z 


^pj 


V /=i+l 


p,j'’ 


(4) 


where p = odd/even with .S'odd,; = ‘^ 1 , + S^yenj - 

S\. + The order parameters Dodd and Deven are each 

non-zero in one phase and vanish in the other. More specif¬ 
ically, Dodd vanishes in the rung-singlet phase and Devon van¬ 
ishes in the Haldane phase. Also, Ocd vanishes in the rung- 
singlet and Haldane phases. The behaviour of the three order 
parameters can thus also be used as a means of locating the 
quantum critical point and the range of each phase. 


yj.^y±/yx, yx^^lyx- (5) 

Consequently if the CD phase were to exist in the region 
0.36 < 7j_ < 0.4 for Jx = 0.2, then there must also be phase 
boundaries in the corresponding region 1.8 < 7j. < 2.0 for 
7x = 5. Yet no further phase boundaries have been observed 
in this region, which is consistent with the evidence for a sin¬ 
gle phase boundary, with no intermediate phase, between the 
rung-singlet and Haldane phases. It thus seems plausable to 
rule out the existence of the CD phase in the antiferromagnetic 
region on symmetry arguments alone. 

It appears then that the only contentious point remaining 
for the cross-coupled spin ladder with antiferromagnetic cou¬ 
plings is identifying the order of the transition along the phase 
boundary. It is clear that along this line a singlet-singlet 
gap vanishes between the gapped Haldane and rung-singlet 
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phases. The early analytic and numerical results were consis¬ 
tent with the transition being first order for both weak and 
strong interchain couplings . However, later work sug¬ 
gested the transition was continuous for weak interchain cou¬ 
pling, becoming first order for stronger couplingJ^ii^ It has 
also been concluded that the transition is first order along the 
whole linei^ Here we have seen that the transition is continu¬ 
ous at 7x = 0.2, = 0.383. This supports there being some 

crossover point(s) on the phase boundary separating lines of 
continuous and first order transitions. Such a point is known 
to exist for the zig-zag ladder model, or equivalently the J\-J 2 


spin chain, which exhibits a single line of first order transitions 
for weak couplings crossing over to second order transitions 
for strong couplings.— The precise clarification of this issue 
for the cross-coupled spin ladder deserves further study. 
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